home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dlaqge.z / dlaqge
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAAQQQQGGGGEEEE((((3333FFFF))))                                                          DDDDLLLLAAAAQQQQGGGGEEEE((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLAQGE - equilibrate a general M by N matrix A using the row and scaling
  10.      factors in the vectors R and C
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DLAQGE( M, N, A, LDA, R, C, ROWCND, COLCND, AMAX, EQUED )
  14.  
  15.          CHARACTER      EQUED
  16.  
  17.          INTEGER        LDA, M, N
  18.  
  19.          DOUBLE         PRECISION AMAX, COLCND, ROWCND
  20.  
  21.          DOUBLE         PRECISION A( LDA, * ), C( * ), R( * )
  22.  
  23. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  24.      DLAQGE equilibrates a general M by N matrix A using the row and scaling
  25.      factors in the vectors R and C.
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      M       (input) INTEGER
  30.              The number of rows of the matrix A.  M >= 0.
  31.  
  32.      N       (input) INTEGER
  33.              The number of columns of the matrix A.  N >= 0.
  34.  
  35.      A       (input/output) DOUBLE PRECISION array, dimension (LDA,N)
  36.              On entry, the M by N matrix A.  On exit, the equilibrated matrix.
  37.              See EQUED for the form of the equilibrated matrix.
  38.  
  39.      LDA     (input) INTEGER
  40.              The leading dimension of the array A.  LDA >= max(M,1).
  41.  
  42.      R       (input) DOUBLE PRECISION array, dimension (M)
  43.              The row scale factors for A.
  44.  
  45.      C       (input) DOUBLE PRECISION array, dimension (N)
  46.              The column scale factors for A.
  47.  
  48.      ROWCND  (input) DOUBLE PRECISION
  49.              Ratio of the smallest R(i) to the largest R(i).
  50.  
  51.      COLCND  (input) DOUBLE PRECISION
  52.              Ratio of the smallest C(i) to the largest C(i).
  53.  
  54.      AMAX    (input) DOUBLE PRECISION
  55.              Absolute value of largest matrix entry.
  56.  
  57.      EQUED   (output) CHARACTER*1
  58.              Specifies the form of equilibration that was done.  = 'N':  No
  59.              equilibration
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAAQQQQGGGGEEEE((((3333FFFF))))                                                          DDDDLLLLAAAAQQQQGGGGEEEE((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              = 'R':  Row equilibration, i.e., A has been premultiplied by
  75.              diag(R).  = 'C':  Column equilibration, i.e., A has been
  76.              postmultiplied by diag(C).  = 'B':  Both row and column
  77.              equilibration, i.e., A has been replaced by diag(R) * A *
  78.              diag(C).
  79.  
  80. PPPPAAAARRRRAAAAMMMMEEEETTTTEEEERRRRSSSS
  81.      THRESH is a threshold value used to decide if row or column scaling
  82.      should be done based on the ratio of the row or column scaling factors.
  83.      If ROWCND < THRESH, row scaling is done, and if COLCND < THRESH, column
  84.      scaling is done.
  85.  
  86.      LARGE and SMALL are threshold values used to decide if row scaling should
  87.      be done based on the absolute size of the largest matrix element.  If
  88.      AMAX > LARGE or AMAX < SMALL, row scaling is done.
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.